Quantum oscillations in topological superconductor candidate Cuo.25Bi2Se3 
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Quantum oscillations are generally studied to resolve the electronic structure of topological insula- 
tors. In Cuo.25Bi2Se3, the prime candidate of topological superconductors, quantum oscillations are 
still not observed in magnetotransport measurement. However, using torque magnetometry, quan- 
tum oscillations (the de Hass - van Alphen effect) were observed in Cuo.25Bi2Se3 . The doping of Cu 
in Bi2Se3 increases the carrier density and the effective mass without increasing the scattering rate 
or decreasing the mean free path. In addition, the Fermi velocity remains the same in Cuo.25Bi2Se3 
as that in Bi2Se3. Our results imply that the insertion of Cu does not change the band structure of 
BiaSeg. 

PACS numbers: 71.18.+y, 74.25. Ha, 74.25. Jb 



Topological insulators and topological superconductors 
are new families of materials with novel electronic states 
[iHS]. Bi2Se3 is a topological insulator (TI) that has 
attracted special interest since Cu-intercalation between 
quintuple layers of Bi2Se3 can induce superconductivity 
below 3.8 K [7]. Unconventional superconductivity aris- 
ing from these topological materials may support Ma- 
jorana fermions and provide a platform for topological 
quantum computation |10H12j . It is proposed that this 
electron-doped topological insulator Cua;Bi2Se3 is a topo- 
logical superconductor (TSC) [13^ . Instead of a bulk insu- 
lating gap and a surface Dirac state, a TSC has a full pair- 
ing gap in the bulk and a topologically protected gapless 
Andreev bound state consisting of Majorana fermions on 
the surface. Specific-heat measurement later elucidated 
that Cuj.Bi2Se3 is a bulk strong-coupling superconduc- 
tor with a full pairing gap |14j . Furthermore, from a 
point-contact spectroscopy study, the Ando group has 
discovered that this Cua;Bi2Se3 exhibits a surface An- 
dreev bound state [H]. All these analyses of possible 
superconducting gap functions have allowed the conclu- 
sion that CuxBi2Se3 is indeed a TSC. This is a big step 
towards the goal of identifying the new topological phase 
of matter. However, the observation of quantum oscil- 
lations, which are generally studied to resolve the elec- 
tronic structure of topological materials, is still missing 
in Cu:rBi2Se3, the leading candidate of TSC. For a TSC, 
as well as the TIs such as Bi2Se3 and Bi2Te3, it is neces- 
sary to observe quantum oscillations in order to resolve 
Landau level quantization and yield direct measurement 
of the bulk and surface electronic states [TBlEn] . As a re- 
sult, an exact measurement of the effective mass and the 
scattering rate of this TSC material is still controversial 

[HISIHll]. 

We solved the problem by measuring the magnetic 
torque of a single crystalline Cuo.25Bi2Se3. Quantum os- 
cillations measured in torque magnetometry, known as 
the de Hass - van Alphen effect, are observed in the high 
field state of the fully superconducting Cuo.25Bi2Se3. 
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FIG. 1: (color online) Quantum oscillation observed by 
torque magnetometry. Magnetic torque r is plotted as a func- 
tion of l/fioH. The lower left panel shows the sketch of the 
measurement setup, where the magnetic field is applied to the 
sample with a tilt angle (j> relative to the crystalline c axis. 
The Fast Fourier Transform (FFT) plot of the torque signal is 
shown in the upper right panel. To generate a FFT plot from 
the T— 1/ noH data, the polynomial background is subtracted 
before the FFT. 

Compared with the oscillation pattern in Bi2Se3, the 
Fermi surface of Cuq. 25812863 becomes larger and stays 
ellipsoidal. We resolve the effective mass and the scat- 
tering rate in the crystal ab plane. The effective mass 
of Cuo. 25612863 is found to be enhanced to 0.19 rUe (free 
electron mass), compared with 0.14 rrif. from the bulk 
state of Bi2Se3. Moreover, the scattering rate and the 
Fermi velocity stay the same in the Cuq. 25612863, as those 
in Bi2Se3. 

Cuo. 25612863 single crystals were grown by slow cooling 
of melted stoichiometric mixtures of high purity elements 
Bl (99.999 %), Cu (99.99 %) and 8e (99.999 %) from 
850° C to 620 °C in a sealed evacuated quartz tube. The 
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crystals were then quenched at 620 °C in cold water, 
resulting in easily cleaved crystals with shiny mirror like 
surfaces. However, the silvery surfaces turn golden after 
one day of exposure to air. The chemical formula of 
Cuo. 25612863 was determined according to the mole ratio 
of the starting elements used in the crystal growth. 

Torque magnetometry was applied to measure the 
anisotropy of magnetic susceptibilities of samples under 
external magnetic fields pm26j . With the tilted mag- 
netic field H confined to the x-z plane, and M in the same 
plane, the torque \s t = M x H = {M^H^ - M^H^)y. 
In the normal state of samples, where the magnetization 
is strictly linear in H, it is convenient to express the ob- 
served torque as a function of the susceptibility difference 
as follows. 



(a) 6 



( b ) , 



= {Xz - Xx)H^H^ ^ Axi?^ sin cos ( 



(1) 



where </> is the tilt angle of H away from z, and 
= Xz — Xx is the magnetic susceptibility anisotropy. 
In this letter, the crystalline c axis of the single crystal 
Cuq. 25612863 is parallel to the z axis. 

In our experimental setup, a Cuq. 25612863 single crys- 
tal is glued to the tip of a thin metal cantilever, as shown 
in the lower panel of Fig. [T] The magnetic torque t is 
measured capacitively. Torque r is plotted as a func- 
tion of the inverse of the applied magnetic field fioH, 
displayed in Fig. [T] at temperature T = 300 mK and 
(f) ^ 11°. Oscillatory pattern in r is periodic in l/figH, 
reflecting the quantization of the landau levels. For met- 
als, the oscillating frequency Fs is determined by the 
cross section area A of the Fermi surface, as following. 



^5 - T^A. 



(2) 



In the upper right panel of Fig. [T] the Fast Fourier Trans- 
formation (FFT) of the r shows a single peak in the am- 
plitude spectrum with a corresponding Fs ~ 325 T. 

The dHvA effect is also observed in the normal state of 
our Cuq. 25612863 samples. One may suspect an extreme 
picture that phase separation within the samples forms 
two different domains: a Cu-rich region which is super- 
conducting and another Cu-deficient region where simply 
the undoped 612863 produces a quantum oscillation sig- 
nal. The first challenge is to figure out whether the oscil- 
latory pattern comes from the Cu-doped superconducting 
phase or from the undoped 612863 phase. To answer this 
question, we first present the evidence of zero-resistance 
and the Meissner effect that are consistent with earlier 
reports from other groups. Then, we compare the dHvA 
effect of Cuo. 25612863 and 612863. We found that the Cu 
doping changes the quantum oscillation pattern, by in- 
creasing the carrier density as well as the effective mass. 




FIG. 2: (color online) The temperature T dependence of the 
sample resistance R. (Panel a) and magnetic susceptibility 
X (Panel b). (Panel a)Zero resistance is observed below 1.2 
K. (Panel b) Volume magnetic susceptibility x is measured at 
field = 2 G in the crystal ab plane, in both the Zero-Field- 
Cooled (ZFC) and the Field-cooled (FC) conditions. From the 
ZFC curve, the nominal superconducting fraction is found to 
be around 35 %, consistent with the fraction found in earlier 
reports [I4ll22] . 



but it does not change the Fermi velocity or the scatter- 
ing rate. We conclude that the quantum oscillation is 
indeed coming from the Cu-doped phase. 

The four-probe resistance R of the sample is shown as 
a function of T in Figj2]ja). As T decreases, R starts 
to drop at 3.3 K, decreases quickly at T = 3 K, and 
becomes zero below T — 1.2 K. The low field magnetic 
susceptibility was also measured with a Quantum Design 
MPM8 magnetometer. Figj2jb) shows the volume mag- 
netic susceptibility x = M/^iqII at 1.8 K < T < 6 K 
aX H — 2 G in-plane magnetic field. In the Field-cooled 
(FC) run, the sample was cooled down to 1.8 K under 
the 2 G external field and x measured during the 
warm up. 6y contrast, in the Zero-Field-cooled (ZFC) 
run, the sample was cooled down at zero external field. 
In both runs, rapid decrease of x was observed at T ~ 
3.5 K. Moreover, at the lowest achieved T ~ 1.8 K in 
the MPMS system, x ~ -0.35 suggested a 35% super- 
conducting volume, which is comparable to other reports 
[m [22]. This considerably large superconducting frac- 
tion suggests that our Cuq. 25612863 sample is indeed high 
quality single phased crystal. Moreover, the observation 
of the zero resistance implies that a large scale phase 
separation does not occur in our single crystals. 

Further evidence shows the effect of Cu doping 
by comparing the dHvA effect of the superconductor 
Cuo. 25612863 and the topological insulator 612863. Fig. 
|3] shows the oscillating period Fg at low tilt angles mea- 
sured in Cuo. 25612863 (Panel a) and 612863 (Panel b). 
First of all, Fs for Cuq. 25612863 is much larger than that 
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of the undoped sample, which suggests the Cu doping 
indeed adds carriers into the electronic state. For both 
samples, the oscillation from the bulk state shows only 
one oscillation frequency, suggesting a single ellipsoidal 
Fermi pocket in both samples. Compared with Bi2Se3, 
Cuo.25Bi2Se3 shows the following three features. 
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FIG. 3: (color online) Angular dependence of the oscillat- 
ing period Fs is compared between (Panel a) Cuo.25Bi2Se3 
and (Panel b) Bi2Se3, Solid lines are the fits based on a sin- 
gle ellipsoidal Fermi surface. An example of the oscillating 
magnetic torque is shown in the inset of Panel b for Bi2Se3. 

1. Larger Fermi pocket. The angular dependence of 
Fs provides an estimate of the Fermi pocket size. The 
Fg — 4> curve n Cuq. 25612863 is consistent with an ellip- 
soidal Fermi surface with the cross section A^z about 
4.02 nm^^. Base on these results and Eq. [2J we obtain 
kf, = k^p = 0.97 nm-\ and fc|, = 1.3 nm^^. Thus, the 
bulk carrier density n = ^^kpkpkp is 4.3 x 10^^ cm~'^. 
This carrier density is similar to the estimate based on 
the Hall effect fl4! and ARPES [21]. For the control sam- 
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pie Bi2Se3, the sinusoidal fitting gives kp 
nm^^ and kp = 1.2 nm~^. The bulk carrier density n is 
calculated to be 1.8x 10^^ cm~'^. This number is larger 
than those generally seen in results based on quantum 
oscillations [18], which suggests that there are quite a 
number of carriers caused by defects. 

We note that not only the carrier density is increased in 
the Cu doped topological superconductors, but also the 
electronic state is closer to a three dimensional sphere. 
The aspect ratio kp/kp of Bi2Se3 is about 1.6, consis- 
tent with ref. [18]. By contrast, kp/kp is about 1.3 
in Cuq. 25612863. The doping of Cu changes the relative 
ratio of Fermi ellipsoid and makes it closer to a sphere. 

2. Slightly increased effective mass. We now focus on 
the temperature dependence and magnetic field depen- 
dence of the oscillating amplitude of magnetic torque as 
the magnetic field is along the c-axis. The oscillating 
torque dr is obtained by subtracting a monotonic back- 
ground from the raw torque data. In metals, the first 
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FIG. 4: (color onhne) Temperature dependence of oscillating 
magnetic torque. (Panel a) Magnetic torque after subtract- 
ing a polynomial background dr is plotted as a function of 
1/H with ~ 4° at selected T between 0.3 K and 30 K. 
Curves at different T are displaced for clarity. (Panel b) In 
Cuo.25Bi2Se3, temperature dependence of an oscillating am- 
plitude dr yield the effective mass m = 0.194me. The dingle 
plot is shown in the inset, giving the Dingle temperature To 
= 23.5 K. (Panel c) Comparison plot in Bi2Se3. The T de- 
pendence of dr at H = 16.5 T yields m = 0.140me, and the 
Dingle plot at 300 mK shows To = 23.9 K. 



harmonic of the oscillating magnetic torque is well de- 
scribed based on the LK formula [27] , and dr is propor- 
tional to the thermal damping factor Rt and the Dingle 
damping factor Rd, as follows. 



Rt = aTm*/Bsmh{aTm,*/B) 
Rd — exp(— aTom/B) 



(3) 



where the effective mass m — m*me, the Dingle temper- 
ature Td — h/2'KkBTs. Ts is the scattering rate , rrig 
is the free electron mass, and a = 27^^fcB"^e/e^ ^14.69 
T/K [H]. 

Fig. Qa) displays the oscillating amplitude dr vs. 
l/HoH at selected T between 300 mK and 30 K. Off- 
sets are added to these curves for clarity. The tilt angle 
(/) ~ 4° . All the following fitting results are for H closely 
along the c -axis. 

First, under a fixed field, the T dependence of normal- 
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ized dr is determined by the effective mass m*me- Fig. 
^h) shows the dr vs. T at if = 16.07 T. The oscihating 
ampHtude dr is normahzed by drg, the amphtude at the 
lowest T = 300 mK. Using the thermal damping Rx for- 
mula in Eq. [3j a fit shown in a line yields the effective 
mass m = 0.194 mg. 

By comparison, the same analysis using the oscillating 
torque dr at H = 16.5 T in Bi2Se3 is shown in Fig. |4](c). 
The Rt fitting yields m = 0.140 mg, which agrees with 
earlier reports [TB]. The comparison shows that the ef- 
fective mass increases gently in the topological supercon- 
ductor candidate Cuq. 25612863 . However, our result does 
not agree with the large mass-enhancement suggested by 
heat capacity measurements |14j . 

3. Similar Fermi velocity and scattering rate. Besides 
adding more carriers and increasing the effective mass, 
the Cu doping does not change the band structure of 
the parent compound Bi2Se3. This can be shown as we 
calculate the Fermi velocity vp — fikp/'ni. Based on the 
effective mass values and fc|,, we obtain vp = 5.8 x 10^ 
m/s for Cuo. 25612863, and vp — 5.7 x 10^ m/s for Bi28e3. 
8ince vp determines the slope of the E — k dispersion at 
the Fermi surface, the relatively similar values of vp in 
both samples suggested the added electrons by Cu doping 
enters the same mobile bulk band in Bi2Se3. 

Further analysis of H dependence of the oscillat- 
ing amplitude gives a similar scattering rate in both 
Cuo.25Bi28e3 and Bi28e3. The inset of Fig. ^h) dis- 
plays the "Dingle plot" at 300 mK in Cuo.25Bi2Se3. The 
fit yields the Dingle temperature Tp, = 23.5 K, which 
gives the scattering rate rg — 5.2 x 10""'^^ s. Thus, the 
mean free path I = vpTg = 30 nm. In Bi2Se3, the same 
fitting shown as the Dingle plot in Fig. |4]jc) yields Tp) — 
23.9 K, Ts = 5.1 X 10"" s, and I = 29 nm. 

TABLE I: Parameters in Cuo. 25612863 and Bi2Se3 



Cuo. 25812863 Bi28e3 

Fs 325 T 150 T 

n 4.3x 10^'' cm"^ 1.8x 10^'' cm"^ 

k%- 0.97 nm-i 0.69 nm"^ 

fe|,/fe^ 1.3 1.6 

m 0.194 me 0.140 

Vp 5.8 x 10® m/s 5.7 x 10** m/s 

TS 5.2 X 10"^" s 5.1 X 10"" s 

I 30 nm 29 nm 



Discussion 

Table |T] summarizes the analysis results of the dHvA 
effect in Cuo.25Bi28e3 and Bi28e3. The doping of Cu in- 
creases the carrier density and the effective mass, but 
does not affect several key band structure parameters es- 
pecially Vp. As Vp determines the slope of the energy 
dispersion at the chemical potential, the unchanged vp 
implies that the added carriers by Cu doping go into the 



same conductive band in the undoped 612863. Note that 
kp increases by more then 40 % in a Cuq. 25612863, which 
should leads to a large increase in vp for a quadratic 
band. Our observed unchanged vp suggests a rather lin- 
ear Dirac-like band. 

The ratio of To to Tc also reflects the relationship be- 
tween the scattering ratio A — l/rg — 2TrkBTp)/h and 
the superconducting gap. Our result shows that To = 
23.5 K is much larger than Tc — 3.5 K. Assuming the 
6C8 gap 2 A — 3.5A;bTc, we find the dimensionless ra- 
tio /cbTd/A ^ 3.8. Therefore, the superconductivity in 
Cuq. 25612863 occurs in the dirty limit. Further, if the 
A/T(7 ratio follows the 6C8 theory, the Pippard length 
= Hvp/ttA is 2.3 /im, which is about two orders of 
magnitude larger than the mean free path /. As a re- 
sult, the coherence length is mainly determined by the 
mean free path and the impurity scattering. We note 
that the ab— plane coherence length determined by the 
upper critical field is 13.9 nm |14j . similar to the mean 
free path. 

Topological odd-parity superconductivity is proposed 
[T^ to exist in Cuq. 25612863. Experimental supports in- 
clude a recent observed zero-bias conductance peak in 
the point-contact spectrum [15]. For a general pairing 
based on the 6C8 theory, regular s-wave pairing is pro- 
tected from the nonmagnetic impurities whereas 
odd-parity pairing is easily destroyed by the impurity 
scattering [551 [SO]- Our observation of the unchanged 
vp implies that the basic band parameter remains the 
same as Cu is doped to the topological insulator Bi28e3. 
Reference [JT] suggests the strong spin-orbit locking in 
this band structure makes the odd-parity pairing robust 
against impurity scattering. 

In conclusion, quantum oscillation is observed in the 
magnetization of Cuo. 25612863. Based on the quan- 
tum oscillating pattern, only one Fermi pocket exists in 
Cuo. 25612863. The doping of Cu to the topological insu- 
lator Bi28e3 increases the carrier density and the effective 
mass without damaging the scattering rate and the mean 
free path. The Fermi velocity stays the same after the 
Cu doping, which implies the band structure is not af- 
fected by the insertion of Cu. Therefore, the analysis of 
the quantum oscillations has shown another important 
piece of evidence that Cuq. 25612863 is likely a topological 
superconductor. Further studies of quantum oscillations 
in other Cu doping levels will reveal how the Cu dop- 
ing tunes the band structure and the fate of topological 
superconductivity. 
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